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On the Poisson Approximation to Photon Distribution for Faint Lasers 
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It is proved, that for a certain kind of input distribution, the strongly binomially attenuated photon number 
distribution can well be approximated by a Poisson distribution. This explains why we can adopt poissonian 
distribution as the photon number statistics for faint lasers. The error of such an approximation is quantitatively 
estimated. Numerical tests are carried out, which coincide with our theoretical estimations. This work lays a sound 
mathematical foundation for the well-known intuitive idea which has been widely used in quantum cryptography. 



1. Introduction 

The security of Quantum Key Distribution 
(QKD) is based on the non-cloning principle of 
an unknown quantum state pQ. In the implement 
of QKD based on BB84 protocol one expect 
that each pulse contains only one photon. If not, 
the eavesdropper can acquire information using 
beamsplitter attack [3] without exposing its exis- 
tence. However, since an ideal single photon state 
is difficult to prepare, practically, faint laser pulse 
with ultra-low mean photon number is used as 
a convenient realization of pseudo-single photon 
source 0]. 

By letting a laser source pass through a strong 
attenuator we get faint laser pulse. For security 
concern the mean photon number in each faint 
laser pulse is kept very small (about 0.1). But 
still, there is a small probability of having more 
than one photon in each pulse. A precise estima- 
tion of this unwelcome probability is crucial for 
security analysis in QKD systems (2). 

In the literature, the photon number in faint 
laser is treated as Poissonian distributed. It is 
all right if the input laser before attenuation is 
Poisson. However, practically we may have in- 
put laser whose photon number statistics is not 
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Poisson If this laser is used as input, the at- 
tenuated laser may not be Poisson either. But 
there is a common belief that no matter what 
distribution the input laser is, if we attenuate it 
into a faint laser with sufficient small mean pho- 
ton number, then Poisson distribution would be 
a good approximation of photon number distri- 
bution in the faint laser pulse. So far, however, 
this claim has not been mathematical rigorously 
proved, which is the motivation of this work. 

We proved, that for a certain kind of probabil- 
ity distribution, after the binomial decay trans- 
formation, which is a mathematical description 
of laser attenuation jS], the decayed distribution 
can well be approximated by a Poisson distribu- 
tion provided that the expectation of the decayed 
distribution is sufficient small. It gives a theoret- 
ical validation of the above claim, i.e., generally 
we can use a Poisson distribution to approximate 
the photon distribution in faint laser and the er- 
ror of the approximation could be neglected. 

2. Preliminary 

Consider N independent particles (photons in 
laser pulses) passing through an attenuator. Each 
particle has a probability of 77 (0 < i] < 1) to 
penetrate the attenuator. We define X to be the 
number of particles before decay (Input); and X v 
to be the number of particles after decay (Out- 
put). X and X n are random variables taking val- 
ues in the natural number system N, and their 
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probability mass functions (PMF) are P(N) and 
P». 

We can connect P(n) and Pr]{n) via the bino- 
mial decay transformation jS] 

00 /AA 

P »= E Jl B (l-/^(jV). (1) 

It is easy to check that 
P v (n) > 0, Vn G N, 
and 



X») = £ P(iV) = 1. 



71=0 



N=0 



So is indeed a PMF. 

The binomial decay transformation establishes 
a relation between photon number distribution 
before and after the attenuation. The remaining 
of this section gives some properties of it. 

Lemma 1. Let E(X) and E{X rj ) denote the ex- 
pectation of the discrete random variables X and 
Xri, respectively. Then E{X r} ) = i]E(X). 

Proof. By the definition of expectation and 
Eq. (JTJ, 

OO 

E(X V ) = X> P » 

71=0 

OO OO , . 

- E"E ( n )v n (i-v) N - n p(N) 

71=0 N=n \ ' 

°° N /AA 

00 

= VNP(N) = r]E(X). 
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Define P A (n) the Poisson PMF with parameter A, 
i.e., P x {n) = e~ A A n /n!, n G N j7E|. We will show 
that the binomial decay transformation preserves 
the Poisson character. 

Lemma 2. Suppose P(N) = P tl (N), then we 
have P n (n) = P^(n). 



Proof. From Eq. JlJ, 



PvH = E (^)v n (l-v) N - n P"(N) 



E 



n\ ^ (N — n)\ 

N=n v ; 



[m(i - v)] 



N- 
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3. Poisson Approximation 



□ 



From Lemma |21 and Lemma ^ we know that 
if P(N), the photon number distribution in the 
laser before attenuation, is Poisson with param- 
eter //, then the decayed distribution, P ?7 (n), is 
also a Poisson whose parameter is ry/i, with 77 be- 
ing the attenuating coefficient. However, prac- 
tically P(N) may not be a Poisson distribution 
[5] If so, P v (n) would not be a Poisson distri- 
bution. Nevertheless, in QKD the faint laser is 
treated as Poisson distributed. The reason of do- 
ing this is based on the common belief that any 
input distribution would reduce to Poisson dis- 
tribution provided that the attenuation is strong 
enough. Next we justify it quantitatively. 

Proposition. In the case of faint laser, the de- 
cayed distribution can be approximated as a Pois- 
son distribution, 

P n {n) PS P A H, neN, 

where 

X = E(X V ) = rjE(X) < 1. 

More concretely, if we expand P A (n) into Taylor 
series of X: 

P A (0) = l-A + ^-0(A 3 ), (2a) 

P A (1) = A-A 2 + 0(A 3 ), (2b) 

P A (2) = y-0(A 3 ), (2c) 

P x (n) = 0(A 3 ),n>3. (2d) 
Then we have 

P„(0) = 1 - A + y + C(X)X 2 - D {X)X 3 , (3a) 
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P I) (l) = A-A 2 -2C(X)A 2 + J Di(X)A 3 , (3b) 



P v (2) = — +C{X)\ 2 -D 2 (X)\\ (3c) 



]T Pv( n ) = i D o(X) + D 2 (X) - D^X)] A 3 , (3d) 

n— 3 

where 

Var(X) - E(X) 



C(X) = 
and 



2E(xy 



< Di{X) < D(X) = ^±JL, (i = 0, 1, 2). 

M(X) = E [X(X - l)(X - 2)] is the 3rd facto- 
rial moment of X. 

Proof. The generating function of X is 

oo 

G(z)= ^?(JV)/,z£l. 

N=0 

Taking the n-th order derivatives of G(z) with 
respect to z yields, 



G {n) {z) = N(N-l) ■ ■ ■ {N-n+\)P{N)i 



N- 



N=n 

Let z = 1, one has, 
XI 



G (n) (l) = E 



{X-n)\_ 

For n = 0, 1, 2, 3 we have 

G(l) = 1, 

G'(l) = E(X), 

G"(l) = Var(X) + [E(X)} 2 -E(X), 

G"'(l) = £[X(X-l)(X-2)] = M(X). 

From Eq. 0J, 

p » = E (^"(i-i)^^^) 

N—n \ U ' 



V N(N-l) ■ ■ ■ (N-n+l)(l-r 1 ) N - n P(N) 



N=n 



I]" 



G^(l-r]). 



Expanding G^(l — 77) into Taylor serials, in the 
case of n = 0, one has, 

P„(0) = G(l)-r 1 G , (l)+^-G"(l)-^G"'(l-9 r ] ) 



1-A- 



A 2 , 7ar(X) - x2 G"'(l-M,3 



2 j b(a:) ;; 



6£(a:) 3 



= 1 - A + — + G(AQA 2 - A)(*)A 3 , 
where 8q G [0, 1] and 



< D (X) = 



G"'(l - 9 Q v) ^ M(X) 



6E(Xf 



< 



E(xy 



= D(X). 



Analogously, Eq. 13 bp and Eq. (|3cj) can be de- 
rived. Finally, apply the relation 

00 

E*n(») = 1 - [^(°) + ^(1) + p v(n > 



n— 3 



one yields Eq. (|3d|) . 



□ 



Remark 1. If C(X) and D(X) is not too big, 
A 2 G(Af ) and A 3 D(X) can be ignored when A < 1. 
By comparing Eqs. (|2*aJl - lj2cHl with Eqs. (jBliJl - ipkl)! 

we can see that P v (n) is well approximated by the 
Poisson distribution P x (n). 

Remark 2 . The approximation of Poisson distri- 
bution is an asymptotic result in the limit A —> 0. 
If A is large, this approximation will be broken, 
as can be found from the numerical Example ^ in 
the following. 

Remark 3. The approximation error A(n) = 
P n {n) — P x (n) can be written as 



A(0) = A 2 G(AT) + 0(A 3 ), 

A(l) = -2A 2 G(AT) + 0(A 3 ), 

A(2) = A 2 G(AT) + 0(A 3 ), 

A(n) = 0(A 3 ), n>3. 



(4a) 
(4b) 
(4c) 
(4d) 



Here C(X) = [Var(X) - E{X)]/[2E{X) 2 } is de- 
termined by P(N), the input distribution only. If 
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Var(X) = E(X), then C(X) = and the error 
decreases to (3(A 3 ). On the other hand, for some 
singular input distribution, C(X) is so big that 
Pfjin) can no longer be approximated by Pois- 
son. In Example|3|we give a typical example that 
Poisson approximation fails. 

In QKD we use faint laser to simulate the single 
photon source. For security analysis, it is impor- 
tant to estimate P v (n > l\n > 0), the probability 
that a pulse contains more than one photon 0|. 
According to our estimation, 



P v (n > l\n > 0) = 



1-P„(0)-P„(1) 
1-P„(0) 

f + C(X)X 2 



A- -C(A)A 2 



Here we have neglected the A 3 and higher order 
terms. We further simplify it by removing the A 2 
terms in the denominator, which gives 



P n (n >l\n> 0) 



\+C { X) 



A. 



(5) 



If the input distribution is Poisson, then 
C(X) = 0. After we attenuate it to faint laser 
that contains an average of 0.1 photon in each 
pulse, which means A = 0.1, we would have 
P v (n > l\n > 0) = 0.05, i.e., each pulse has about 
5% chance to contain more than one photon. 

From Eq. J5J we can see that, the risk of a 
QKD system rises as C{X) grows. In Example 
13 we use an ill-shaped input distribution whose 
C(X) = 9.11. After we attenuate it to A = 0.1, 
then P v (n > l\n > 0) « 1. So for this big a C(X), 
the QKD system would be totally unsecured since 
almost every pulse contains at least two photons. 

On the other hand, if the input distribution 
satisfies Var(X) < E(X) then C(X) < 0. One 
can expects the attenuated faint laser be more se- 
cure than a Poisson laser because we get a smaller 
P v (n > l\n > 0). The last row of table^gives an 
example of negative C(X). However, the possible 
smallest C{X) is — ^E(X)^ 1 . For practical input 
laser, E(X) is so big that — ^E{X)~ 1 can be ig- 
nored. This means, for attenuated faint laser, a 
Poisson distribution is almost as good as one can 
expects. 
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Figure 1. If attenuation is not strong enough, the 
decayed distribution may not be approximated by 
Poisson distribution. In this example the input 
random variable takes integer values from to 
1000, r\ is set to 0.1, and the expectation after 
decay is A = 48.85. 



4. Numerical Examples 

Following we give some numerical simulations 
and the results coincide with our theoretical esti- 
mations quite well. 

Example 1 . First we show that if the mean value 
of the decayed distribution A = T]E(X) is not suf- 
ficiently small, then in general, approximation us- 
ing Poisson distribution fails. Wc choose a ran- 
dom variable X whose PMF takes the shape in 
Figure^and set rj = 0.1. In this case A « 48.85, 
and it can be observed that P v (n) is far away from 
Poisson distribution. 

Example 2. Using the same input random vari- 
able X as Example ^ but we take r/ = 0.001, 
which means that our new attenuator is 100 
times stronger than the old one. Now we have 
A = 0.4885 and the binomial decayed PMF P n (n) 
is close to a Poisson distribution, as indicated in 
Figure |2 We further decrease rj to 0.0002, then 
A = 0.0977 and Figure [3] shows a perfect match 
between the decayed distribution and Poisson dis- 
tribution, which claims strong support for the 



5 




Figure 2. Using the same input distribution as 
Example ^ but we replace the attenuator with a 
100 times stronger one with r\ = 0.001 and A = 
0.4885. Only the probability value of the first few 
n's are plotted since the others are almost zero. 



common belief that photon number in faint laser 
pulse can be treated as Poisson distributed. 

Example 3. For some ill-shaped input distribu- 
tions, C(X) could be very large. If this happens, 
even if A is small, Poisson approximation could 
still fail. As an example, we construct a singular 
input distribution of which P(X = 1) = 0.95 and 
P(X = 1001) = 0.05. In this case, C{X) = 9.11. 
If we attenuate it to A = 0.1, as Figure 01 shows, 
P x (n) fails to converge to P v {n). 

To quantitatively test the approximation error, 
Eqs. H4a jl -I|4d fl . we choose different input distribu- 
tions and computer their attenuated distribution. 
We adjust the attenuating coefficient t] to keep 
A = 0.1. The numerical result, which is listed in 
Tabled supports our theoretical estimation very 
well. 
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Figure 3. Like Example we use a even stronger 
attenuator with r] = 0.0002 and A = 0.0977. 
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Table 1 

Numerical validation of the approximate error, Eqs. (|4all - 14dll . It can be observed that A(0) « A(2) ~ 
X 2 C(X), which goes on well with Eq. |ga} and Eq. (|4H1) : And A(l) ps -2A 2 C(A), which also agrees with 
Eq. gcjl. 



X 2 C{X) 


A(0) 


A(l) 


A(2) 


A(3) 


A(4) 


0.0045 


0.0039 


-0.0073 


0.0028 


0.0005 


0.0000 


0.0030 


0.0027 


-0.0050 


0.0020 


0.0003 


0.0000 


0.0018 


0.0017 


-0.0032 


0.0013 


0.0002 


0.0000 


0.0011 


0.0010 


-0.0019 


0.0008 


0.0001 


0.0000 


0.0005 


0.0005 


-0.0010 


0.0004 


0.0000 


0.0000 


-0.00016 


-0.00014 


0.00027 


-0.00011 


-0.00001 


-0.00000 




Figure 4. A typical example when error is large 
even if A is small. The input distribution is con- 
structed as P(X = 1) = 0.95 and P(X = 1001) = 
0.05. We decay it to A = 0.1 and it turns out that 
P x {n) fails to converge to P v {n). 



